Archaeobotanical evidence supports indigenous cucurbit long-term use in the Mesoamerican Neotropics

The squash family (Cucurbitaceae) contains some of the most important crops cultivated worldwide and has played an important ecological, economic, and cultural role for millennia. In the American tropics, squashes were among the first cultivated crop species, but little is known about how their domestication unfolded. Here, we employ direct radiocarbon dating and morphological analyses of desiccated cucurbit seeds, rinds, and stems from El Gigante Rockshelter in Honduras to reconstruct human practices of selection and cultivation of Lagenaria siceraria, Cucurbita pepo, and Cucurbita moschata. Direct radiocarbon dating indicates that humans started using Lagenaria and wild Cucurbita starting ~ 10,950 calendar years before present (cal B.P.), primarily as watertight vessels and possibly as cooking and drinking containers. A rind directly dated to 11,150–10,765 cal B.P. represents the oldest known bottle gourd in the Americas. Domesticated C. moschata subsequently appeared ~ 4035 cal B.P., followed by domesticated C. pepo ~ 2190 cal B.P. associated with increasing evidence for their use as food crops. Multivariate statistical analysis of seed size and shape show that the archaeological C. pepo assemblage exhibits significant variability, representing at least three varieties: one similar to present-day zucchini, another like present-day vegetable marrow, and a native cultivar without modern analogs. Our archaeobotanical data supports the hypothesis that Indigenous cucurbit use started in the Early Holocene, and that agricultural complexity during the Late Holocene involved selective breeding that encouraged crop diversification.


El Gigante Rockshelter
El Gigante (14° 13ʹ N 88° 03ʹ W, 1300 m above sea level) is a large rock shelter located along the western ridge of the Estanzuela River in the highlands of western Honduras (Fig. 1).The site was excavated in 2000 and 2001 24,25 .The occupation at the site spans great time depth with deposits dating from the Terminal Pleistocene through the Postclassic/Colonial periods.Bayesian chronological modeling of radiocarbon dates from the rock shelter suggest 16 phases of occupation 26 .The initial occupation took place during the Early and Middle Esperanza phases between 10,985-10,705 and 10,705-10,210 cal B.P., respectively, and was characterized by the presence of stemmed bifaces and avocado, squash and maguey remains.After a gap of 190 years, the Late Esperanza (10,020-9,520 cal B.P.) featured wide fluted projectile points as well as ciruela (Spondias), squash, maguey, and avocado remains 27 .The subsequent phase of occupation, Marcala, was episodic.The Early Archaic (Early Marcala 1: 8945-8520 cal B.P., Early Marcala 2: 8090-7865 cal B.P.) and Middle Archaic (Middle Marcala 1: 7570-7265 cal B.P., Middle Marcala 2: 7140-6960 cal B.P., and Middle Marcala 3: 6665-6365 cal B.P.) habitations are characterized by the abundance of ciruela, avocado, squash, and maguey remains.A long occupation hiatus took place between 6365 and 4440 cal B.P. The Late Marcala phase (4400-4025 cal B.P.) corresponding to the regional Late Archaic is characterized by the appearance of maize (Zea mays), which is found along with tree fruits, including ciruela and avocado as well as a low abundance of squash.Following a second hiatus of 530 years, the Early and Middle Formative occupation (Early Estanzuela phase: 3500-3185 cal B.P.; Middle Estanzuela: 2985-2470 cal B.P.) features the appearance of ceramics.The Late Formative occupation (Late Estanzuela phase 1: 2225-2080 cal B.P.; Late Estanzuela phase 2: 1980-1810 cal B.P.) contains large quantities of maize and the appearance of beans (Phaseolus spp.).During the Classic period (Early Classic: 1690-1475 cal B.P., Late Classic: 1295-865 cal B.P.), maize continues to constitute a main component of the assemblage, and avocados, squash, and wild fruits increase in abundance.The occupation of the rock shelter lasted through the Postclassic/Colonial eras (545-310 cal B.P.), periods when the site was rarely used.During the most recent occupations, acorns were common, but maize and gourds were rare in the assemblage.

Direct radiocarbon dating
Direct accelerated mass spectrometry (AMS) radiocarbon ( 14 C) dating of desiccated cucurbit rinds, stems, and seeds throughout the 11,000-year occupation of El Gigante has revealed a long history of cucurbit use by humans (Table 1).The earliest analyzed Cucurbita rinds date from the Esperanza phase and were directly AMS 14

Peduncles
A total of 16 Cucurbita peduncles was recovered, of which seven were identified as C. pepo and seven as C. moschata (following 14, 15, 49, 50).C. pepo stems exhibited seven to nine major ridges that had a mean basal diameter of 19.7 ± 0.5 mm (Fig. 3c).C. moschata peduncles were characterized by the presence of five major ridges, a pentagonal base (Fig. 3d), and a mean basal diameter of 13.9 ± 1.7 mm.
To better understand the variability in the archaeological assemblage, we compared the archaeological C. pepo seeds from El Gigante to a dataset of modern wild and domesticated C. pepo seeds that we measured (Supplementary Table S1).Principle component analysis (PCA) of morphometric traits allow us to separate modern variants of domesticated C. pepo seeds from those of wild C. pepo (C.pepo var.fraterna, C. pepo subsp.texana and C. pepo var.ozarkana) (Fig. 5a).The first two components collectively explain 78% of the variance in the data (61.5% and 16.8%, sequentially).Principal component 1 discriminated between total seed length and length to width ratio, while principal component 2 separated seeds by maximum width and seed margin width.Measurements of archaeological seeds recovered from El Gigante demonstrate that most cucurbit seeds are comparable in size to modern domesticated C. pepo seeds.Only six of the El Gigante seeds are sufficiently small to fall within the confidence interval ranges of wild C. pepo seeds and specifically, the varieties C. pepo var.ozarkana and C. pepo var.fraterna.
To evaluate if the ancient C. pepo seeds recovered from El Gigante are similar to modern varieties of domesticated and wild C. pepo, we performed a second PCA using the cucurbit cultivar measurement data reported in Ref. 40 .Principal component 1 (63% of total variance) discriminated among seeds by total length and maximum width whereas principal component 2 (35.2% of total variance) separated seeds according to their length to width ratio (Fig. 5b).The biplot shows high variability within the El Gigante archaeological seeds and the fact that the first two components captured 98% of the variance attests to the strength of the relationships among the cases and variables.Ten archaeological seeds were morphometrically similar to modern acorn, scallop and crookneck varieties of wild C. pepo subsp.texana.Sixteen seeds overlap with modern seeds of cocozelle, vegetable marrow and zucchini varieties of C. pepo subsp.pepo.Finally, 17 of the archaeological seeds are not associated with any variety of modern C. pepo in particular.
A multivariate analysis shows that only one archaeological seed of C. moschata from El Gigante is morphometrically similar to its modern equivalent (Fig. 5c).Principal component 1, which accounts for 43.6% of the variance, separates seeds by total length, maximum width, and mass.Principal component 2, accounting for 30.3% of the variance, separates seeds by surface area and length to width ratio.Three of four archaeological seeds fell outside the confidence intervals of modern C. moschata.Archaeological seeds have low principal component 1 scores, which are reflective of small seeds (small length and small width), but high principal component 2 scores, indicating that they have a broader surface area than most of the modern seeds.
A final PCA comparing ancient and modern L. siceraria seeds shows that the first component (62.8% of the variance) differentiates seed surface area, total length, and mass.The second component (26.8% of the variance) is related to maximum width and length to width ratio (Fig. 5d).Most of the El Gigante Lagenaria seeds are distributed along the left side of the first component, showing their smaller size in relation to modern seeds.To assess shape variation, we calculated elliptical Fourier descriptors (EFDs) (Fig. 6).The results largely show that Cucurbita species are clearly discriminated from Lagenaria siceraria on the LDA1-LDA2 biplot (94% of the total variation).Furthermore, the LDA discriminates the seeds of domesticated C. pepo from the wild C. pepo species (C.pepo var.fraterna, C. pepo var.ozarkana, and C. pepo subsp.texana).It is notable, however, that there is significant overlap in seed shape between C. pepo and C. moschata and between C. pepo var.ozarkana and C. moschata.
The geometric morphometric analysis indicates that the shape of the archaeological gourd and squash seeds from El Gigante vary considerably.Of the 21 archaeological seeds morphologically identified as C. pepo, the elliptical Fourier transform (EFT) method correctly classified only three specimens as C. pepo.The remaining archaeological seeds were allocated to the wild C. pepo group (seven as C. pepo var.ozarkana, eight as C. pepo var.texana, and three as C. pepo var.fraterna).Only one of the archaeological C. moschata seeds was correctly classified as C. moschata, while the other two seeds were attributed to the C. pepo group.Finally, most of the archaeological Lagenaria seeds (n = 5) were correctly classified as L. siceraria, but three were assigned as C. pepo and one as C. pepo var.ozarkana.In light of the variability inherent in cucurbit seeds, this method has limited power to add to morphological classification in this case.

Discussion
The chronology, size, and geometric morphometric analyses of Cucurbita and Lagenaria macrobotanical remains from El Gigante Rockshelter reveal intriguing patterns related to human reliance, use, and cultivation of cucurbits for food and technology during the last ~ 11,000 years.Key changes in plant selection and use correspond with broader subsistence and settlement changes spanning the Esperanza and Marcala phases, a time span that represents the long transition from plant foraging to low-level cultivation.
Ancient people inhabiting the area around El Gigante consistently used L. siceraria and Cucurbita through the entire sequence of occupation of the site.Initially, Cucurbita rinds likely gathered from the wild were significantly thinner than rinds recovered from subsequent periods when cultivated varieties were probably utilized.Research has shown that an increase in rind thickness is associated with the presence of domesticated squash in contrast to the brittle rinds of wild squash forms 12 .El Gigante Cucurbita rinds recovered from the earliest occupation levels were < 2 mm thick, suggesting that humans used wild Cucurbita species during this time.In contrast, Cucurbita rind thickness increases during the Late Archaic (2.5 mm, range 0.5-5.1),indicating the introduction of domesticated Cucurbita to the region ~ 4000 years ago.A similar change in rind thickness over time has been documented in the Guilá Naquitz cave, where Cucurbita rinds showed an increase in thickness from 0.84 mm (range 0.5-1.6 mm) to 1.15 mm (range 0.5-2.0mm) between ca.10,500 and 8500 years ago 12 .In contrast, rinds of C. maxima recovered from the Pampa Grande archaeological site in Argentina show high variation in thickness through time (mean 4.09 mm, range 2.08-6.64mm), suggesting that different varieties were cultivated and managed for different purposes (e.g., some as food and others as containers) 41,42 .
Given that a large percentage of Cucurbita rinds from these earlier occupations were burned, it seems likely that mobile hunter-gatherers used some cucurbits as cooking containers.We also suggest that some of these early, thinner cucurbits were roasted directly in proximity to cooking fires for consumption, potentially focusing on the squashes that had a lower concentration of bitter cucurbitacin than the rest of the fruits.In contrast, likely cultivated Cucurbita with more durable and thicker rinds could have been used to store food or water.The recovery of a Cucurbita rind stained with red pigment suggests that squash containers also may have been decorated vessels 42,43 .The recovery of Lagenaria rinds and seeds spanning the occupation of El Gigante Rockshelter provides further evidence that gourds were used consistently during the early Holocene in Mesoamerica.AMS dating of a Lagenaria rind recovered from Unit 18, level 39 corresponding to the early Esperanza phase produced a date of 10,935 cal B.P. (11,150-10,765 cal B.P.).The date is considerably older than the most ancient specimens of bottle gourd previously reported, which were from Guilá Naquitz, Oaxaca (10,021 ± 105 cal B.P.; 10,198-9781 cal year BP) and Little Salt Spring, Florida (10,015 cal B.P.) 8,44 .Lagenaria was likely one of the first field crops in El Gigante and the light-weight fruits were used as containers even after ceramic technology was introduced to the region during the Early Formative period (~ 3600 cal B.P.).
Lagenaria rinds show little change in thickness over time.It is notable that El Gigante Lagenaria rinds exhibit a thickness (1.9 ± 1 mm) that overlaps with rinds of wild African L. siceraria (< 2.0 mm) 3 .El Gigante rinds are also consistently thinner than those recovered from other archaeological contexts in the Americas, including the Ocampo (3.2 mm), the Guilá Naquitz (4.6 mm), and the Coxcatlán (3.8 mm) caves in Mexico, and Ancón (4.8 mm) in Peru 3 .The domestication status of the recovered Lagenaria remains to be assessed in detail due to the absence of morphological characterizations in wild and domesticated Lagenaria species.These results highlight the importance of conducting further comparative studies of modern and archaeological plant remains alongside paleogenomic analyses to assess the domestication syndrome in Lagenaria and its possible center of domestication in south-central Mesoamerica.
Morphometric studies of variation in seed size improve our understanding of how domestication shaped plant diversity 45 .Refined methods of morphometric analysis are critical if we are to fully understand the processes involved in the domestication and spread of these cultigens.Geometric morphometric analysis is a complementary approach that can be used to assess changes in seed shape that may be associated with domestication.This method has been useful in the analysis of several crop species, including olives, grapes, barley, and date palms [46][47][48][49] .For instance, seed morphology has been used to distinguish between pips of wild and domesticated grape (Vitis vinifera) varieties in Italy 50 .Our results indicate that domesticated C. pepo show discernable differentiation in seed size and shape from wild species (C.pepo var.fraterna, C. pepo subsp.texana, and C. pepo var.ozarkana).Furthermore, distinct groups of C. pepo and C. pepo visualized with PCA and EFDs confirm existing species distinctions based on morphological data.We argue that the high level of phenotypic variability in cultivated C. pepo is a consequence of varietal diversification and reflects the long-term history of cultivation associated with the selection of specific morphological traits, such as fruit size and shape.These morphological differences may be explained by divergent selection pressures.For instance, wild C. pepo species were subject to environmental constraints that seemed to have minimized seed size, while domesticated C. pepo were subject to ongoing human selection related to cultivation practices that shaped particular phenotypes.
In contrast to fruit size, it is likely that seed shape was not a direct target of human selection in C. pepo.The evolutionary role of seed shape in cucurbits is not well understood, but shape changes could be related to a pleiotropic effect that determines seed persistence in the soil, seed germination and dispersal 51 .Some empirical studies have shown that long-term persistent species tend to have spherical seeds, while transient species produce flattened and/or elongated seeds 52,53 .Genetic drift associated with the restricted gene pool of breeding can also have essentially random implications for non-selected traits such as seed shape.
The shape analysis was not able to differentiate seeds of C. moschata from subspecies of C. pepo.Modern seeds of C. moschata showed a high overlap with both the domesticated C. pepo and the wild C. pepo var.ozarkana.The results expose the limitations of using shape analysis for identifying C. moschata seeds and highlights the importance of including other qualitative characters for identification.
Multivariate morphometric analyses show that most of the C. pepo seeds from El Gigante fall within the size parameters defined for modern types of domesticated C. pepo seeds.Paris and Nerson 39 have also shown a relationship between seed size and fruit size in C. pepo.Our analysis suggests that some of the archaeological seeds are similar to domesticated present-day varieties, including the cocozelle and zucchini varieties, which are characterized by long and uniformly cylindrical fruits, and the vegetable marrow variety, which is distinguished by short round fruits with a slightly broad apex 39 .
The seeds of C. pepo that are similar in size to wild C. pepo subspecies could have been the product of genetic exchange between crops and wild cucurbit species in the area.Gene flow and crop-to-wild hybridization have been reported to be common in C. pepo due to pollen transfer between plants by bees 40,54 , and hybrid plants have been shown to produce highly variable fruit sizes 54 .
In contrast to the morphometric results, shape analysis via linear discriminant analysis (LDA) allocated most of the El Gigante C. pepo seeds to wild forms.These results could suggest the presence of an unknown autochthonous variety or a separate subspecies that has since become extinct, for instance the as-yet-unknown wild progenitor of C. pepo ssp.pepo cultivated types.They also suggest that the region in which El Gigante is located may have been an important setting of squash diversification.A molecular analysis could help to better address squash diversity via aDNA analysis.However, the difficulty of differentiating wild from domesticated seeds in archaeological contexts could be associated with the preservation of the seed remains.Even though this analysis only included complete and well-preserved seeds, it is possible that the desiccation process might have altered the shape of some seeds.

Conclusions
The archaeobotanical record from El Gigante Rockshelter suggests early indigenous foragers utilized wild cucurbits for generations before the introduction of cultivated varieties.Identifications of maize pollen and microcharcoal particles recovered from wetlands in Lago Yojoa ~ 5400 cal B.P. and in the Copan Valley ~ 4300 cal B.P. suggests that the adoption of farming in Mesoamerica was a relatively early and protracted process 55,56 .In correspondence with these studies, domesticated butternut squash seems to have been initially introduced in www.nature.com/scientificreports/El Gigante during the Late Marcala phase, between 4400 and 4025 cal B.P., yet the utilization of wild cucurbits and other tropical plants persisted over time.In fact, as indigenous farmers began to substantially modify their surrounding tropical ecosystems, they were able to successfully incorporate domesticated crops such as maize and squash while still relying on and likely managing other resources such as tree crops (e.g., Persea americana, Pouteria spp., and Spondias spp.) for their subsistence 26 .The significance of El Gigante lies in both the length of its chronology and its extensive well-preserved plant assemblage.The rock shelter thus has great potential to both broaden and reshape our understanding of longterm human adaptations to domestication and early agriculture.The analytical techniques we used (digital morphometrics, multivariate statistics, and outline analysis based on EFDs) proved to be powerful tools for assessing variation in rind thickness and seed size and shape.Indeed, geometric morphometric approaches combined with traditional biometric methods of length and width measurements have allowed us to distinguish among cucurbit species, specifically domesticated C. pepo, wild C. pepo, and Lagenaria.We note, however, that the classification accuracy of C. moschata from C. pepo based on EFDs is more limited compared to qualitative morphological attributes.
The study of the plant remains recovered from El Gigante Rockshelter substantially contributes to improving our understanding of Indigenous human-plant coevolution and the dispersal of key cultigens 57,58 .As part of broader research, this paper presents an important piece of the tropical Americas domestication puzzle.

Recovery of macrobotanical remains
El Gigante is one of six dry rock shelters in Mesoamerica that contain well-preserved macrobotanical assemblages that include both desiccated and carbonized remains 57 .The plant remains analyzed in this study were either hand-collected in situ during excavations or recovered through sieving (mesh size 4.5 mm and 500 μm) sediment from Unit 18. Well-preserved cucurbit (Cucurbita and Lagenaria) remains in the assemblage include rinds, peduncles, and seeds.Cucurbita rinds (n = 448) were distinguished from Lagenaria rinds (n = 387) based on diagnostic cross-section cell structure (Supplementary Fig. S1).Cucurbita rinds exhibit isodiametric cells with a regular configuration, a feature that enabled us to distinguish them from Lagenaria rinds 12 .Cucurbita peduncles were identified at the species level based on the shape of the peduncle.Wild Cucurbita presents relatively small peduncles (< 13 mm in diameter) with circular outline; domesticated C. pepo peduncles are relatively bigger (> 14 mm in diameter) and are characterized by the presence of alternating ten major bridges with an angular or pentagonal outline at the fruit attachment, while C. moschata peduncle presents five ridges, with regular furrows and flaring at the base 12,13,59,60 .

Direct radiocarbon dating
To further verify the age of the cucurbit assemblage, we conducted direct AMS 14 C dating of 28 samples, including some of the earliest specimens in the assemblage.Samples were pretreated using conventional acid-base-acid methods and subsequently combusted, graphitized, and measured at the accelerated mass spectrometer of the Pennsylvania State University Radiocarbon Laboratory (PSUAMS).In addition, direct 14 C dating of 310 additional (non-cucurbit) plant specimens allowed us to confidently assign the majority of cucurbit specimens to clarify chronological phases 26 .

Modern seed material
To assess differences in seed shape and size, we relied on modern seeds from the USDA Germplasm Resources Information Network (GRIN).Samples included accessions from wild Cucurbita species (C.pepo subsp.texana  S4).We also relied on previously published measurements of additional varieties of domesticated subspecies and cultivar groups of C. pepo 39 .

Seed size and shape analysis
Digital images of the modern and archaeological seed samples were acquired using a digital camera attached to an Olympus SZX16 microscope with a 300-dpi resolution.Modern and archaeological seeds were photographed dorsally.Seeds were not photographed laterally because most of the archaeological seed coats were empty, resulting in lateral shape distortion.Digital images of the seeds were processed and measured using ImageJ v. 1.42 61 .The allometric characteristics of maximum length, maximum width, border thickness, and area were recorded in millimeters.Seeds were also weighed using a digital scale with a measuring error of 0.01 g.The full measurement datasets are included in Supplementary Tables S1 and S2.
Individual digital pictures were converted to binary images (i.e., seed image on a black background) in ImageJ.Each binary image was converted into chain code along the perimeter of each seed to create a harmonic series using SHAPE v. 1.3 62 .Seed shape was quantified using outline analysis based on EFDs to assess for differences in overall seed shape within species.Chain code contours were posteriorly converted to normalized EFDs for Fourier analysis.

Data analysis
To assess differences in rind thickness between phases of occupation, we performed a one-way ANOVA for Lagenaria and Cucurbita, respectively.Differences in seed size between modern and archaeological seeds of C. pepo, C. pepo, C. moschata, and L. siceraria were assessed using a Principal Components Analysis (PCA).A second PCA was performed to compare similarities in seed size between seeds of different varieties of C. pepo from Paris and Nerson 39 and archaeological seeds from El Gigante.A Ward's agglomerative hierarchical classification using Euclidean distance was used to assess the morphological relationship of C. pepo seeds from El Gigante to other seeds recovered from paleontological and archaeological sites from North America (Supplementary Table S3).The cluster analysis was conducted using maximum length, maximum width, and length/width ratio.A Linear Discriminant Analysis (LDA) was performed to explore differences in seed shape between modern seeds of C. pepo, C. pepo, C. moschata, and L. siceraria.LDA is a statistical classification method used to evaluate shape variation of seeds by maximizing differences between predefined groups (i.e., species, subspecies and varieties) compared to intragroup variation.
Figure 1.(a) Location of El Gigante Rockshelter and archeological sites with Cucurbita measured seeds used in the analysis.(b) Percentages of macrobotanical remains of Cucurbita and Lagenaria and other field crops from El Gigante Rockshelter, Honduras.

Figure 2 .
Figure 2. Variation of rind thickness of (a) Cucurbita and (b) Lagenaria recovered across occupational phases from El Gigante Rockshelter, Honduras.Different letters indicate significant differences (P < 0.05).The sample size is indicated below each boxplot.

Table 1 .
Direct radiocarbon dates of desiccated cucurbits samples from El Gigante Rockshelter, Honduras.